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Abstract — Performance characteristics of orifice plates for 
variation in different geometrical parameters are analyzed 
with the help of computational fluid dynamics (CFD) for non- 
standard conditions with water as working fluid. In the 
present work, A CFD tool, ANSYS, FLUENT has been used to 
predict the variation of Coefficient of Discharge (C d ) and the 
analysis for standard sharp edge orifice meter has been done 
for different plate thicknesses (3mm, 5mm, 10mm and 15mm) 
in a pipe of 50mm diameter and the effect of Pipe diameter on 
Coefficient of Discharge (C d ) has been studied in detail. The 
value of Coefficient of Discharge depends on the type of flow, 
pressure tappings, contour of the obstruction and it is a 
function of Reynolds number. The results and outcomes of the 
study are presented with velocity and pressure contours. The 
effect of non-standard conditions on Coefficient of Discharge 
(Cd) has been discussed in detail. 

Index Terms — Coefficient of Discharge (Cd), CFD, Reynolds 
number, ANSYS, FLUENT, Orifice Meter 

I. Introduction 


To scientifically analyze them is a challenge. This attempt 
is to see if CFD tool can be used for accurately calculating 
Cd in such cases and to develop a methodology that could 
be used for solving CFD problems under non-standard 
conditions. The results are positive if the CFD tool is 
properly handled. 

Extensive work has been done in previous studies 
and many have analyzed the effects of various parameters 
on the Coefficient of Discharge. Singh et.al [l]conducted a 
detailed study on the effects of small diameter orifice 
meters, with different beta ratios and plate thickness on Cd 
using GAMBIT and ANSYS Rahman et.al [4] have done an 
experimental work in hydraulic lab for different beta ratios 
using a setup and they analyzed and gave a relationship of 
beta ratio with Re and Cd. Naveenji Arun et.al [6] have 
studied the flow of non-Newtonian fluids using CFD and 
plotted the variation of Cd with different concentrations of 
fluids and obtained that C d increases exponentially with 
increase in Reynolds Number. 


Flow measurement is essential part in many industries 
today. The flow rate is to be found out at a check point 
which can be achieved by the device called a flow 
measuring device or a flow meter. The most common type 
used is the Differential Pressure type flow meters, 
sometimes also called the Head-Loss or Obstruction type 
flow meters. In any Obstruction type flow meter, the flow 
rate is calculated by measuring the pressure drop over an 
obstruction inserted in the path of the flow. 

An orifice plate is a thin plate with a hole in it 
which is placed inside the pipeline for the determination of 
flow rate of the flowing fluid. The orifice plate may also be 
used for reducing pressure or for restricting a flow in case 
of restriction orifice plate. When a fluid passes through the 
orifice; its pressure builds up slightly upstream of the pipe 
and the fluid is forced to pass through the hole. The velocity 
increases and the fluid pressure decreases as the area 
reduces due to the presence of obstruction. A little 
downstream of the orifice the flow reaches its point of 
maximum convergence called the vena Contracta where the 
velocity reaches its maximum value. This is where the 
pressures are noted down by the tappings. 

We know that not always standard dimensions and 
procedures can be used for an application. Industries have 
processes which demand or require non-standard operating 
range. There may be aging factors due to years of usage and 
erosion, operating under high pressures, non-accurate 
manufacturing or requirement of non-standard operating 
range for the industrial applications. 
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II. WORKING PRINCIPLE OF ORIFICE METER 
ASSEMBLY: 


The working principle of any obstruction type 
flow meter is based on Bernoulli’s Equation given by 

Pi v l2 p 2 V 2 2 

f Z 1 + 7T" f z 2 + TT- 

Pg 2g pg 2g 

The volumetric flow rate can be expressed by the 
equation: 
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Where dland d2 are the diameters of the pipeline 
and the orifice opening, pi and p2 are upstream and 
downstream pressures, then the flow rate can be obtained 
using above equation by measuring the pressure difference 
(pl-p2). 

The flow expression obtained in the above 
equation is not accurate in a practical case since it does not 
take into account the real fluid effects and hence an 
empirical factor called discharge co-efficient (Cd) is 
incorporated. Thus the equation becomes 
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III. DISCHARGE COEFFICIENT OF ORIFICE PFATE: 

Cd is defined as the ratio of the actual flow to the 
theoretical flow and is always less than one. There are 
many reasons due to which the actual flow rate is less than 
the ideal one. Mainly the assumptions of frictionless flow 
and incompressible flows are not always valid. The extent 
of friction depends on the Reynolds number (Re). As a 
result, the correction factor Cd<l, has to be incorporated. 
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Calculation of Cd can be done using the formula 
specified in ISO-5167. The discharge coefficient, Cd is 
given by the Reader-Harri s/Gallagher Equation which is 
purely empirical and is based on extensive experimental 
data. Their correlation is given by, 

c d = 0.5961 + 0.0261 p 2 - O.2160 8 
+ 0.00052 l[10 6 p/Re D ] 

+ [0.0188 + 0.0063A] (3 35 [10 5 /Re D ] 

+ [0.043 + 0.080e _1 ° Ll 

T 0 4 

— 0.123e -7Ll ][l — 0.11A] — 1 —— 


— 0.031[M 2 — 0.8M 11 2 ]p 1 - 3 
If D<71.12mm (2.8 in), the following term should 
be added to the equation. 
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Where D is the diameter of pipe in mm 
Diameter ratio, p =d/D 
Re D =Pipe Reynolds number 
A= [19000p/Re D ] 0.8 
M2=2L2/(l-p) 

LI and L2are lengths from the faces to tappings and depend 
on the type of taps used. 

The above equation has its own limitations. The 
equation is completely empirical and is evaluated only by 
experimental results. 

It is only applicable when pipe diameter is within 
50mm<D< 1000mm when diameter ratio is within 0.2< 
p<0.75, plate thickness between 0.005D and 0.02D, bevel 
angle of 45°+ 15°. Further the upstream edge of the orifice 
hole should be sharp and a radius of curvature not greater 
than 0.0004d. 

The geometry and dimensions should be as per 
ISO 5167 standards for the equation to be applicable. For 
any other non-standard cases, the equation is invalid. 


IV. CFD & Mathematical Model 

CFD is the most popular tool that is used 
nowadays to analyze a problem. CFD gains its popularity 
due to its capability to produce good results without 
conducting any experiments. 

The governing equations used by the CFD code Fluent 
that has been used in the present study are 

• The continuity equation or the conservation of 
Mass Equation for two dimensional 
incompressible flows. 

dp d(pii) d(pv) d(pw) _ 

dt dx dy dz 

• The Navier Stokes Equation or the conservation of 
Momentum Equation for two dimensional 
incompressible flows. 

du dp /d 2 u d 2 u d 2 u\ 

P dt = PS_ + ^ + 3+ + ) 

dv dp /d 2 v d 2 v d 2 v\ 

P dt = P8_ a^ +P Vdx i + dy I+ dz I J 

In the case of turbulent flow, additional equations 
need to be solved depending on the turbulence model that is 
used in the computational analysis. The details of these 
models are not explained in detail since they are available 


in the standard literature. 

The standard k-s model (Launder and Spalding, 
1974) is used in the simulations carried out and the model 
has two model equations, one for k (turbulent kinetic 
energy) and one for 8 (rate of viscous dissipation). 

The standard k-s model uses the following 
governing equations for k and s shown below. 
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Where the eddy viscosity, |i t is given by the following 
equation 

k 2 

Pt = pc^ — 

Where, is a dimensionless constant. 

The k-s model is the simplest turbulence model 
which only requires initial and boundary conditions well 
established and most widely validated turbulence model. It 
provides excellent results for many industry oriented flows. 
But the model produces poor results in case of some 
unconfined flows with large strains, swirling flows, 
recirculation and rotating flows etc. thus the k-co model 
(Wilcox, 1988, 1993a, b, 1994)can be used as the best 
alternative which uses turbulence frequency cd as second 
variable. 


The flow analysis is done using the help of 
commercially available CFD tool, ANSYS, FLUENT 
which uses finite volume approach to discretize the 
problem. 


V. Validation 

Validation is an important step in any problem 
solving technique. In CFD, validation needs to be done to 
prove that the mesh density used, type of mesh used, 
boundary conditions, turbulence model used in the 
simulation are proper and the results obtained are accurate. 
The validation may be done by comparing the results 
obtained by CFD with the experimental results obtained by 
conducting experiments or by comparing it with the 
available standard values in the ISO-5167 and BS-1042 
standards. 

The validation for the present problem was done 
for the flow through standard orifice plate of pipe diameter 
50mm and Reynolds number of 1, 00,000 for beta ratio of 
0.5 with inlet Velocity of 2m/s and exit gauge pressure of 
OPa. 


Fig 1- 


The geometry of the above problem is as shown in 
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Fig 1: Geometry and Boundary conditions of Standard 
Orifice plate. 



Fig 2: Mesh details of Standard Orifice Plate. 

The modeling was done in the ANSYS 
Workbench 14.0. The model was divided into 3 parts for 
efficient meshing. The upstream length was chosen to be 
15D and the downstream length of 20D was incorporated to 
avoid the interference of the boundary conditions. The Grid 
independence test was done starting from 50000 elements 
to 300000 elements to ensure the minimum mesh density 
required to meet the required accuracy of Cd. It was 
observed that after 140000 the value of Cd remained almost 
constant and hence the meshing was done using 
Quadrilateral elements of 1, 40,000elements as shown in 
Fig 2. Fine meshing was done near the orifice plate due to 
the computational limits and fine wall meshing near the 
wall to capture the variations. The necessary boundary 
conditions at the inlet, outlet and wall were applied (Inlet- 
Velocity inlet, Wall- Wall/No slip condition, Axis-Axis, 
Outlet -Pressure outlet). Problem was considered to be axis 
symmetric. Second order upwind discretization and 
SIMPLE scheme was utilized for solving. Residuals or 
convergence criteria of 10‘ 6 was employed with double 
precision solver. Standard k-s model was used in the 
computations. 

The coefficient of discharge (further referred to as 
C d ) for flange and corner taps were calculated using the 
pressure readings obtained at the respective pressure taps in 
CFD. The results obtained in CFD were validated against 
the ISO -5167 standard values and are shown in Table -1 
and the respective velocity and pressure contours and plots 
are shown in Fig 3. 


Turbulent Flow, Re= 100,000 

Cd 

Flange taps 

Corner Taps 

Value of C d from ISO-5 167 
Standards 

0.6082 

0.6069 

Value of C d obtained by CFD 

0.624 

0.624 


Table 1: Validation results for Standard Sharp Edge Orifice 
Plate. 
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Velocity Contour: Max Value :1.20e+01 m/s; Min Value: 0 

m/s 
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Pressure Contour: Max Value : 7.00e+04 Pa, Min Value: - 
4.86e+04 Pa 

Fig 3: Velocity and Pressure Contours Flow through 
Standard Orifice Plate. 

The above results show that the computed values 
of C d are somewhat higher than the values given in 
standards. Further calculations using different turbulence 
models have shown that the choice of model is also very 
important. The computations made with K-co SST model 
gave the values of C d as 0.6051 for flange taps and 0.6054 
for corner taps. Thus the computed and standard values 
agree with each other within the uncertainty limits of 0.6% 
mentioned in the ISO-5167 standards. Thus the 
computational methodology can be considered as validated. 

VI. Results and Discussions 


EFFECT OF PIPE DIAMETER ON C D 

As per ISO 5167, the minimum pipe diameter for 
which the correlation for C d is applicable is specified as 
5 0mm. The effect of pipe diameter on C d when it becomes 
less than 50mm has been studied in this section. The ranges 
of diameters selected for the study are from 50mm to 10mm 
and the diameter ratio is 0.5. Computations have been done 
for two Reynolds numbers namely 200 (Laminar Regime) 
and 100000 (Turbulent Regime).The obtained CFD results 
for various pipe diameters are tabulated in Table-2, and the 
corresponding velocity and pressure contours for laminar 
and turbulent cases are shown in Fig 4. As the diameter of 
the pipe is decreased, the conclusion that one can draw is 
that C d increases gradually. The variation in flange taps is 
more as compared to the variation observed in Corner taps. 
In flange taps, the pressure tappings are located at a 
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distance of 1 inch (25.4mm) and in corner taps the pressure 
taps are located at the plate corners and are independent of 
pipe diameter. Thus, in the case of flange taps the 
downstream pressure tap in relation to pipe diameter moves 
farther away from the orifice plate as the pipe diameter is 
decreased. Hence the pressure recovery will increase with 
decrease in diameter leading to a lower value of pressure 
differential. Consequently the value of C d would 
correspondingly increase. This has been supported by the 
computations. However this effect is subdued in the case of 
corner taps. Hence it can be concluded that the equation 
given in the standards to calculate C d is not applicable to 
the diameters below 5 Omm. The extent of error will keep on 
increasing with decreasing diameter. 


Serial 

No. 

Diameter 
of the 
pipe , D, 

mm 

C d for Laminar 
Flow, Re =200 

C d for Turbulent 
Flow, Re= 100,000 

Flange 

Taps 

Corner 

Taps 

Flange 

Taps 

Corner 

Taps 

1 

50 

0.720 

0.689 

0.624 

0.624 

2 

40 

0.727 

0.723 

0.627 

0.626 

3 

30 

0.730 

0.727 

0.642 

0.632 

4 

20 

0.750 

0.746 

0.695 

0.653 

5 

10 

0.811 

0.770 

0.823 

0.707 


Table 2: Effect of Pipe Diameter on Coefficient of 
Discharge. 


LAMINAR FLOW 



Pipe Diameter D,mm 
Variation of C d for Laminar Flow 
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3.74e+04 Pa 

Fig 4: Contour plots of Non-Standard Orifice Plate 
(D= 1 Omm) for Laminar flow. 
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TURBULENT FLOW 



Pipe Diameter D, mm 
Variation of C d for Turbulent Flow 
Fig 4: Variations of Cd for Laminar and Turbulent cases for 
different pipe Diameters of Standard Orifice Plate. 


Velocity Contour: Max Value : 1.19e+01m/s ; Min Value: 
Om/s 



Contours of Static Pressure (pascal) 


Pressure Contour: Max Value : 5.39e+04 Pa, Min Value: - 
3.25e+04 Pa 
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Fig 5 : Contour plots of Non-Standard Orifice Plate 
(D=10mm) for Turbulent Flow. 



Velocity Vector: Max Value : 1.19e+01 m/s; Min 
Value: 6.98e-03 m/s 

Fig 6: Vectors plots of Non-Standard Orifice Plate 
(D=10mm) for Turbulent Flow. 
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FLANGE 

TAPS 
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Variation of Cd for Laminar Flow 


EFFECT OF ORIFICE PLATE THICKNESS ON C D 

The maximum allowable thickness of the standard 
orifice plate as per standards for a 50mm diameter pipe is 
3mm. In certain industrial applications, the pipeline 
pressure would be very high and to use an orifice plate to 
measure the flow rate in these cases higher plate 
thicknesses would be required in order to ensure structural 
safety. Thus for design considerations, the thickness of the 
orifice plates are increased. The effect of increased plate 
thickness on Cdis studied using CFD tool. The results given 
below are simulated for a 50mm diameter pipe and Beta 
ratio of 0.5 for both laminar and turbulent cases. 
Computations have been made for six plate thicknesses in 
the range 1mm to 15mm. The results are tabulated in table 3 
for the two types of pressure tappings and at two Reynolds 
numbers. 

The corresponding velocity and pressure contours for 
laminar and turbulent cases for a plate thickness of 5mm are 
given in Fig 8. 

It is observed from the tabulated values that the 
computed values of C d are independent of plate thickness as 
long as it is less than 5mm. In the case of thicker plates, the 
values of discharge coefficient increases with increasing 
thickness and the deviations are fairly significant. This 
shows that the correlation given in ISO 5167 could be used 
for plate thicknesses up to 5mm in a 50mm pipeline. 


Serial 

No. 

Thickness 

of the 
orifice 
plate, t, 

mm 

C d for Laminar 
Flow, Re =200 

C d for Turbulent 
Flow, 

Re= 100,000 

Flang 

e 

T one 

Corner 

Taps 

Flange 

Taps 

Corner 

Taps 

1 

1 

0.699 

1 

0.693 

0.620 

0.620 

2 

2 

0.713 

0.711 

0.624 

0.624 

3 

3 

0.720 

0.689 

0.624 

0.624 

4 

5 

0.72 

0.717 

0.624 

0.624 

5 

10 

0.762 

0.758 

0.68 

0.68 

6 

15 

0.777 

0.772 

0.706 

0.700 


Table3: Effect of plate thickness on Coefficient of 
discharge. 


TURBULENT FLOW 



Pipe Diameter D, mm 
Variation of C d for Turbulent Flow 


Fig 7: Variations of Cd for Laminar and Turbulent cases for 
different Orifice Plate Thickness of Standard Orifice Plate. 
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Contours of Static Pressure (pascal) Mar 31. 2015 
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Pressure Contour: Max Value : 6.2/e+04 Pa, Min Value: - 
3.29e+04 Pa 

Fig 8: Contour plots of Non-Standard Orifice Plate 
(t=5mm) for Laminar flow. 
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Pressure Contour: Max Value : 5.84e+04 Pa, Min Value: - 
2.43e+04 Pa 

Fig9: Contour plots of Non-Standard Orifice Plate (t=5mm) 
for Turbulent flow. 
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Fig 10: Vectors plots of Non-Standard Orifice Plate 
(t=5mm) for Turbulent Flow. 


increases. Similarly, if the orifice plate thickness is larger 
than 5 mm (in a pipe of 50mm diameter) the value of Cd can 
be substantially higher. Thus the study has demonstrated 
the utility of CFD in predicting the performance 
characteristics of the orifice plate assemblies. 
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VII. Concluding Remarks : 

A CFD Methodology for analyzing flow through 
an orifice plate assembly has been presented. This 
methodology has been validated by comparing the results 
with the correlations available in the standards for a 
standard concentric square edged orifice plate. As long as 
turbulence model is carefully chosen, accurate estimation of 
C d can be obtained from CFD analysis. 

Further it has been demonstrated that CFD can 
also be used to predict the values of C d even under non- 
standard conditions. It has been shown that if the pipe 
diameter is decreased below 50mm the value of C d 
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